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1 Introduction 



Effective Lagrangian of the dual magnetic Monopole Nambu Jona Lasinio model. 

In our recent publications [1-4] we have shown that the dual Monopole Nambu-Jona-Lasinio 
model (MNJL) with dual Dirac strings is a perfect continuum space— time analogy of Com- 
pact Quantum Electrodynamics (CQED) [5]. As has been shown in Ref.[5] CQED possesses 
the same non-perturbative phenomena as low-energy QCD. The MNJL model is based on a 
Lagrangian, invariant under magnetic £7(1) symmetry, with massless magnetic monopoles self- 
coupled through a local four-monopole interaction [1,2]: 

C{x) = xix^d^xix) + G[x(x)x(x)} 2 - Gi[x{x)-f^x{x)][x{x)^x{x)l (1.1) 

where x( x ) is a massless magnetic monopole field, G and G\ are positive phenomenological 
constants. Below we will show that we have to choose G\ = G/4 for the self consistency of the 
theory in the one loop approximation [3,4]. 

The magnetic monopole condensation accompanies the creation of XX collective excitations 
with the quantum numbers of a scalar Higgs meson field p and a dual-vector field C^. 

For the derivation of an effective Lagrangian the p and C M fields are introduced as cyclic 
variables. 

C(x) = xi^vfdyxix) - V(x), (1.2) 

where V(x) is defined 

- V(x) = x(x)(-grC^x) -k P {x))x{x) - ^p\x) + ^-C M (x)C"*(x). (1.3) 

Now we can show that the vacuum expectation value (v.e.v) of the p field does not vanish. For 
this aim we have to derive the equation of motion of the p field by varying the Lagrangian 
Eq. (|l.l|) with respect to the p field: 

-AJ. = - K X(x)x(x)-—p{x) = 0. (1.4) 

This leads to 

2G 

p(x) = x(x)x(x)- (1.5) 

K 

Taking the v.e.v. of both sides of Eq.( |1.5| ) we get 

2G 2G 

(p(x)) = (x(x)x(x)) = < X(0)x(0)>, (1-6) 

where (x(O)x(O)) is the magnetic monopole condensate. Thus, the non-zero value of the v.e.v. 
of the p field is related to the monopole condensation. In order to deal with a physical scalar 
field, the a-field, we should follow the standard procedure and subtract <p(x)>. This gives 
<j{x) = p(x) — <p(x)>. The v.e.v. of the p-field can be expressed in terms of the mass of 
the magnetic monopole M in the superconducting phase [1—4], <p(x)>= M/k, where M is 
proportional to <x(0)x(0)> [1~4] 

M = -2G(x(0) X (0)). (1.7) 

This is the gap-equation testifying the appearance of massive magnetic monopoles in the super- 
conducting phase, where <x(0)x(0)>7^ 0- As has been shown in Refs.[l,2] it also leads to the 
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suppression of direct transitions between the physical scalar field a and the non-perturbative 
vacuum. 

In terms of the cr-field the Lagrangian Eq,( |1.2[) reads 



C{x) = x{x){i-fd,> - M) X (x) - V(x), 



where now V(x) reads 

- V(x) = xWi-Q-fC^x) - na{x)) X {x) - ^p 2 (x) + ^-C M (x)C"*(x). 
Integrating out the magnetic monopole fields we arrive at the effective Lagrangian 



(1.9) 



A „(.'•) = 4s - ^p\x) + A- c,u ■)( ■''(.»•) 



(1.10) 



with >C(x) e fT defined as 



^eff(^) 



l ( X 



Det(i<9 — M + <3?) 



Det(ifl - M) 



(1.11) 



Here we have denoted $ = —g^C^ — no, and a = p — M/k. 

The effective Lagrangian C c s(x) can be represented by an infinite series 



n \ 



n=l 



1 



M 



n=l 



(1.12) 



The index L means the evaluation of the trace over the Lorentz indices. The effective Lagrangian 
(x) is given by 



"- 1 *xe# 
(2tt)« 



n 



d 4 xed 4 ki _ i ^ . Xl _ , . . _ i k n . x 



1 1 . 
n 167T 2 ' 



7T 2 i 



xtr L 



1 



k M - k 
x ...*(s n _i) 



*(zi) 



1 



M - k - h 
1 



$(x 2 ) 



M 



A: 



■$(x) 



n-l 



(1.13) 



at A;i + A)2 + . . . + fc n = 0. The r.h.s. of Eq.( 1.13 ) describes the one-massive-monopole loop 
diagram with n-vertices. The monopole-loop diagrams with two vertices (n = 2) determine 
the kinetic term of the cr-field and give the contribution to the kinetic term of the C^-field, 
while the diagrams with (n > 3) describe the vertices of interactions of the a and C M fields. In 
accordance with the prescription given in [1,2] the effective Lagrangian C e s{x) should be defined 
by the set of divergent one-massive-monopole-loop diagrams with n = 1,2,3 and 4 vertices. 
The evaluation of these diagrams gives 



Asff (x) 



1 K 



2 8vr 2 
1 

+ - 



J 2 (,M)d tl a{x)& t a[x) - M 



^ + ^(M)-4M^, 2 (M) 



a 2 (x) 



3 



■2M«4 J 2 (M)a 3 (x) - -k 2 ^J 2 {M)(j A {x) 
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48 7T 2 

2 



8vr 2 

J 2 (M)C^(x)C(x)^ 



8vr 2 



+ 
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4G 



16 vr 2 



[Ji(M) + M 2 J 2 (M)] 



C M (x) C^x) 



(1.14) 



where we have defined C^(x) = W(x) - d v C»(x). Then, J X (M) and J 2 (M) are the following 
quadratically and logarithmically divergent integrals 



k 2 ) 



A 2 - M 2 in + 



r ,*,S f d4fi 1 



k 2 f 



ln\l + 



A 2 \ 


A 2 


M 2 j 


M 2 + A 2 




i A 2 


M 2 


1 M 2 + A 2 



(1.15) 



In order to get correct factors of the a and field kinetic terms we have to set [1,2] 



12 7T 



J 2 (M) 



i 7T 



J 2 (M) 



1. 



(1.16) 



So, the coupling constants are connected by the relation n 2 = 2g 2 /3 [1,2]. 

The effective Lagrangian Eq.( 1.14 ) contains a term linear in the u-field. This part of the 
effective Lagrangian leads to direct transitions a — ► vacuum. In the case of a physical <r-field 
such transitions should be suppressed. In order to suppress these transitions we have to impose 
the constraint [1-4] 



1 



MM) = 0, 



G 2vr 2 

where J\{M) can be connected with the monopole condensate [1-4] 

1 



(1.17) 



< x(o)x(o) >= 



4tt 2 



AfJi(M). 



(1.18) 



Inserting Eq.( 1.18| ) into Eq. (|1.17|) we arrive at the gap-equation ( |1.7| ) 

The coefficient in front of the last term in Eq.( [L.14 ) defines the mass of the field: 



M, 



C 



2Gi 



> 7T^ 



[MM) + M 2 J 2 (M)]. 



(1.19) 



Piling up the gap-equation (1.7) and the constraint ( 1.16 ) we recast the effective Lagrangian 
( 1.14j ) into the form 

C cS (x) = ~C^(x)C^(x) + ±M 2 C^x)C»(x) + 



+ l -d,o{x)d»a(x)- l -Mla 2 {x) 



1 +K- 



a[x 



M a 



~dC^(x)dC^(x) + ^M 2 C^x)C»{x) + 



<7(x) 9"a(<r) - ^M 2 cr 2 (x) + AntK^)], 



(1.20) 
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where M a = 2M is the mass of the <7-field and £; n t[o"(x)] describes the self-interactions of the 
(T-field 



dutHx)] = - K M a a 3 (x) - ^K 2 a 4 (x). (1.21) 

As has been shown in Ref.[3,4] the gap-equation, in turn, can be derived in the one-monopole 
loop approximation by using only the Lagrangian Eq.( |l.l|) . When comparing these two gap- 
equations we get the relation G\ = G/4 that reduces the number of input parameters [3,4]. 
Wave function of the non perturbative vacuum. The MNJL model as well as the NJL 
model [6] and the BCS theory of superconductivity [7] possesses a non-trivial non-perturbative 
vacuum with a wave function [2] 



io><*>= n 

p,A=±l 



|0>, (1.22) 



where (3p = p/E^ = p/\/p 2 + M 2 is the velocity of a massive monopole with the mass M, 
and a(°"(p, A) (or b(°"(—p, A)) denotes the creation operator of a massless monopole ( or anti- 
monopole) with a momentum p and helicity A; |0) = |0)(°) is the wave function of the perturbative 
vacuum of the non-condensed phase. The wave function |0)( M ) of the non-perturbative vacuum 
is distinctly invariant under magnetic U(l) symmetry. The former implies that in the condensed 
phase the magnetic U(l) symmetry is not broken [2]. This is the main peculiarity of the MNJL 
model with respect to the dual Higgs model with dual Dirac strings [8-10] , where the magnetic 
U(l) symmetry is broken spontaneously in the superconducting phase. 

The wave function of the non-perturbative vacuum |0)( M ) is invariant under party transfor- 
mations V, V\0) {M) = \o)( M \ 

In order to show that the magnetic monopole condensate (x(O)x(O)} nas a distinct meaning 
of the order parameter we, following the BCS theory of superconductivity [7], can introduce two 
operators possessing different properties under parity transformations 

p A=±l 

O- = 2^Y1 b^{-p,X)a^\p,X),VO^ = -G_. (1.23) 

p A=±l 

As has been shown in Ref.[2] the v.e.v of the 0+ operator per unit volume coincides with the 
magnetic monopole condensate (0+) = (x(0)x(0)} : 



p A=±l p p 

d 3 p f d 4 p 1 M 



f d d p [ <rp 1 M , . , 

= - iM Jwr- iM J im i^w^ = -^ Ji{M)=m)xm ' (1 ' 24) 

where V is a normalization volume. In turn, one can show that (O-) = 0. This confirms the 
parity conservation in the MNJL model. 

Dual Dirac strings. Dual Dirac strings are included in the MNJL model in the form of a dual 
electric field strength E IJ,u {x) defined by [1-4] 



5 



where = X^(t, a) represents the position of a point on the world sheet swept by the string. 
The sheet is parameterized by internal coordinates — oo < r < oo and < a < it, so that 
X^(t, 0) = X^q(t) and X^(t, it) = Xq(t) represent the world lines of an anti-quark and a 
quark [1-4,8-10]. 

The effective Lagrangian of the dual Higgs field a(x) and vector C^(x) fields is then defined 
C eS (x) = - A F^{x)F^{x) + ^MlC^x)C»{x) + 



+ l -d,a{x)d^(x) - l -Mla 2 {x) 



1 +K 



<nx} 

Mrr 



\f, u {x)F^{x) + \mIC,{x)C»{x) + 



+-d u a(x)d fl a(x) 



1 



M^a 2 (x) + £ int [a(x)}, 



(1.26) 



where the field strength F^(x) reads [1-4,8-10]: F^(x) = S^(x) - *C^(x) and *C^(x) is 
the dual version of C^(x), *C^(x) = \e^ vafi C afi (x) (e 0123 = 1). 

Quarks and anti quarks. In the MNJL model a quark and an anti-quarks are point-like 
particles with masses m q = rriq = m, electric charges Q q = —Qq = Q, and trajectories Xq(r) 
and X^(t), respectively, attached to the ends of a dual Dirac string. They are described by the 
Lagrangian 



^free quark ^ ^ "tj j 



dr 



l=q,q 



( dX?{r)dX»{T) 
1 dr 



1/2 



d,T 



g, u \ 5^(x-X t (r)). 



(1.27) 



Substituting Eq.(|j~27l) in Eq.(|L2|) we arrive at the total effective Lagrangian of the MNJL 
model with dual Dirac strings, quarks and anti-quarks 



\f, u {x)F^(x) + ±M*C,(x)C»(x) + 



1 



1 



d^a(x) - ^M 2 a 2 (x) 



1 + K 



Mrr 



~F» u (x)F^(x) + ±M 2 c C,(x)C»(x) + 



+~<^ a{x) d^a(x) - ^Ml<j 2 {x) + C irA [a{x)\ 

"^•^free quark 



(1.28) 



This effective Lagrangian should be used for the evaluation of different observables defined in 
terms of vacuum expectation values related to the magnetic monopole Green functions [1-4]. 
Magnetic monopole Green functions. The n-point magnetic monopole Green function can 
be defined as the vacuum expectation value of the time ordered product of the massless magnetic 
monopole densities Refs.[l-4]: 



G(x\,... ,x n ) = (0|T(x(xi)rix(xi) . . . x(x n ) r nX (O)|0)c 



(1.29) 



where Ti{i = l,...,n) are the Dirac matrices. As has been shown in Ref.[ll] the vacuum 
expectation value Eq.( |1.29| ) can be represented in terms of the vacuum expectation values of the 
densities of the massive magnetic monopole fields xm(x) coupled to the fields of the collective 
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excitations a and 



G(xx, ...,x n ) = (0|T(x(a;i)rix(a;i) . . . x(xn)TnX(x n ))\0)conn. = 

= (M) (0|T^XM(»l)riXM(»l) • • -XMiXn^nXuiXn) 

x expi / d A x{-gx M (x)YxM(x)C u (x) - kxm(x)xm(x)ct(x) + C int [a(x)]})\0)i^. 



(1.30) 

Here |0}( M ) is the wave function of the non-perturbative vacuum of the MNJL model in the 
condensed phase and |0) the wave function of the perturbative vacuum of the non-condensed 
phase. 

The self-interactions £i nt [o"(x)] provide cr-field loop contributions and can be dropped out 
in the tree cr-field approximation Refs. [1-4]. The tree cr-field approximation can be justified 
keeping massive magnetic monopoles very heavy, i.e. M S> Mc- This corresponds to the London 
limit M„ = 2M> Mc in the dual Higgs model with dual Dirac strings [8-10]. The inequality 
M a 3> Mc means also that in the MNJL model we deal with Dual Superconductivity of type E 
[12]. In the tree cr-field approximation the r.h.s. of Eq.( |1.3d ) can be recast into the form 



G(xi, ...,x n )=< 0|T(x(»i)r x(O)|0> c 

= (AQ< 0\T (xm{xi)Tixm{xi) . . . XM{x n )^nXM{x n ) 



expz 



J d A x{ - gXM(x)YxM(x)C u (x) - hxm(x)xm(x)ct(x)}^ |0 . (1.31) 



For the subsequent investigation it is convenient to represent the r.h.s. of Eq.( 1.3l| ) in terms of 
the generating functional of the monopole Green functions [1-4] 



n S S 

G(xi,...,x n ) = T] \ T i x -t -, z [v,v] 

JLA 5ri(xi) or] (xi) 



(1.32) 

77=77=0 



where 77(77) are the external sources of the massive monopole (antimonopole) fields, and Z[rj,rj] 
is the generating functional of the monopole Green functions defined by 

Zfafj] = \\ VxmVxmVC^Vo expi J d 4 x [1 (x) F^ (x) 

+~ Ml C,(x) C»{x) + ~ d,a{x) d»o{x) - \ Ml a\x) 
+XM{x){i 'fd lt -M-g'f C^x) - K(j(x)) xm(x) 

+V(x)xm(x) + XM(x)r](x) + £f rec quark (z) • (1-33) 

The normalization factor Z is defined by the condition Z[0, 0] = 1. 

The paper is organized as follows. In Sect. 2 we discuss the gluon condensate in non- 
perturbative QCD both in the QCD sum rules approach and in lattice simulations. In Sect. 3 
we calculate the dual-vector field condensate and compare the obtained result with the gluon 
condensate of non-perturbative QCD. In the Conclusion we discuss the obtained results. 

2 Gluon condensate of non-perturbative QCD 

For the first time, the condensate of the gluon fields 

{■§-2 G lM G°^(0))^0, (2.1) 
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where g s and (0) (a = 1, . . . , 8) are the quark-gluon coupling constant and the gluon field 
strength, has been introduced as a phenomenological parameter characterizing quantitatively 
non-trivial properties of a non-perturbative vacuum of QCD in the QCD sum rules approach 
[13]. The gluon condensate breaks the dilatation invariance at a scale of energy transferred of 
order Arj ~ 4GeV [14]. This gives a signal to the breaking of chiral symmetry which becomes 
broken spontaneously at a scale of energy transferred of order A x ~ 1 GeV [14] . As has been 
shown in Ref.[15] the contribution of the gluon condensate to the condensate of light quarks 
u, d and s makes up 1/3 of the meanvalue of the quark condensate (uu) = (dd) = (ss) = 
-(0.253 GeV) 3 [15]. 

The numerical value of the gluon condensate 



({^G^(0)G a ^(0)) = (0.331 GeV) 4 



(2.2) 



obtained in Ref.[13] by studying the charmonium channel has been obviously underestimated 
[16-18]. The correct value of the gluon condensate increased by a factor of 4 has been obtained 
in Ref.[19]. The meanvalue of the gluon condensate is equal to [19] 



({^G^(0)G a ^(0)) = (0.458 GeV) 4 . 
In the dilute-instanton gas approximation the gluon condensate has been defined as [13] 



(2.3) 



(J^(O)G^(O) 



16 



inst.+anti— inst. 



dp 
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(2.4) 



where d(p) is the instanton density function. For the SU(N) gauge group the instanton density 
function is defined by [20] 



d(p) 



Cie 



-CoN 



(N-l)\(N-2)\ 



8tt 2 
9 2 S (P) 



2N 



^ 2 / 9 2 Ap) 



(2.5) 



The coefficients C\ and C2 are given by [20] 

2e 5 / 6 



Cl 
C 2 



0.466, 



7T^ 



5 „ 11 1 , „ , v 2 inn 

<n2 + (M2lr) + 7)+ E =uw, 

n=0 



(2.6) 



where 7 = 0.5277 ... is Euler's constant. 

In the absence of quark contributions and in the one-loop approximation the running cou- 
pling constant g(p) is determined [20] 



9 2 (P) 



9s (Po) 



1 + u N iM £n ± 



(2.7) 



P0 



where po = 1/Atj and Ajj is the ultra-violet cut-off. Then, the integral over p defining the 
gluon condensate Eq.( |2,4j ) is infrared divergent and the parameter p c = 1/Ar plays the role of 
the infrared cut-off [13]. According to estimates of Ref.[13] the infrared cut-off should of order 
A R ~ 200 MeV. 
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Thus, in the absence of quark contributions and in the one-loop approximation the instanton 
density function d{p) reads [20] 



d(p) 



-CoN 



(N - 1)\(N - 2)! 
For QCD when iV = 3we obtain 

d(p) = 1.513 (pA v 



( P Au) 11N/3 



8n 2 11 



2N 



-87r 2 / ff 2(Au) 



(2.8) 



.11 



8vr 2 



5 2 (Au) 



11 £n(p Ai 



-87r 2 / ff 2 (Au) 



(2.9) 



Substituting Eq. (|2.9[) in Eq^^) one can see that the integral over p is concentrated near the 
infrared cut-off p c = 1/Ar. Therefore, the result of the integration can be given by the expression 
[13,20] 



1.51310" 



16 



A 



inst.+anti— inst. 



R 



8vr 2 



3.458 10 -3 A R 



5 2 (Ar) 

6 



8tt z 



5 2 (Ar) 



e -8^ 2 / 5 2 (A R ; 
-8vr 2 / 5 2 (A R ) 



(2.10) 



According to lattice formulation of QCD the gluon condensate can be represented in the following 
general form [21] 



^(|^^(0)G^(0) 



A 4 A 2 
^40 + W 04 + W 22 + 



A 4 



A 2 



(2.11) 



where W mn are the numerical coefficients and Aqcd enters to the definition of the running 
coupling constant [19,22] 



9s (P) ^ 



1 



47T 



bo 



(2.12) 



in 



pA QC D 



where b = (ll/3)iV [22]. The typical value of A QCD is A QCD = 100 4- 300 MeV [22]. 

The coefficient W40 describes a perturbative contribution to the gluon condensate, whilst the 
coefficients W40 and W22 are fully non-perturbative. The presence of the term proportional to 
Aqqq/Au differs the expression for the gluon condensate obtained within the dilute-instanton 
gas approximation from the gluon condensate calculated on lattice. The appearance of this 
term is related to power corrections [13,21], whereas the expression Eq.( 2.10| ) corresponds to the 
leading order contribution in power expansion [13]. The gluon condensate given by Eqs.( p.l0Tj 
and (PH) we would compare with the condensate of the dual-vector field which we calculate 
below in the MNJL model. 



3 Dual— vector field condensate in the MNJL model 

In the MNJL model we defined the condensate of the dual-vector field by analogy with the 
magnetic monopole condensate [2-4] 



W 



1 

Z 



D X DxDC»Da 



9 



9 



expi J <fz y- C„(z) (□ + Mq) C v {z) + C v {z) d»*£^(z) 
+ ^a{z)(U + Ml)a(z)+C- int [a(z)] 

-9 Xm{z)^vXm{z) C u (z) - kxm{z)xm{z)<j{z) + XM{z)(i^/ud u - M)xm(z) 
where Z is the normalization constant 

Z = JJJJ D X D X DC^Da expi J d 4 z ^-C u (z){U + M 2 c )C v (z) 

+ C v {z) d»*£^(z) + \ *(z) (□ + Ml) a(z) + C- mt [a{z)} 

-9 'Xm{z)1uXm{z) C v \z) - kxm{z)xm{z)<j{z) + XM{z)(i~i v d v - M) xm(z) 



(3.1) 



(3.2) 



Recall that we are working in the London limit, M CT = 2M> Mq- One can show that in this 
limit the contribution of the a field exchanges to the dual-vector field condensate are insignificant 
and can be neglected with respect to the contributions of the dual-vector field exchanges. As a 
result the integral over the a field can be absorbed by the normalization constant Z 



(£^(*)^(*)) = 1 JJjDxDxDC^ 

/' 



^C, u {x)C^{x) 



expi d z 



C v {z) (□ + M 2 C ) C v {z) + C{z) d»*£^(z) -9Xm{z) 1vX m{z) C v {z) 



+XM(z)(ijvd u - M)xm(z) 



(3.3) 



The integration over the dual-vector field we perform around the Abrikosov flux line. For 
this aim we represent the dual-vector field C M as follows [3,4,9,10] 



C u (x)=C u [£(x)]+c v (x). 



(3.4) 



The field C tl [£{x)\ is the Abrikosov flux line induced by a dual Dirac string obeying the equation 
[1-4,8-10] 

p + M*)C v [£(x)] = -&»S^(x) 



(3.5) 
(3.6) 

r) 4 Ml -fc 2 -i0 " (3 " ?) 

The field c M (x) stands for the quantum fluctuations of the dual-vector field around the Abrikosov 
flux line [2-4,9-10]. Denoting F^ v [£{x)] = C^ u [£{x)} and subtracting a trivial contribution 
proportional to F fll/ [£(x)]F^ u [£(x)] we obtain 



possessing the solution 

C^[£(x)\ = - J d 4 x f A(x-x',M c )d»*£^(x f ) 
where A(x — x', Mc) is the Green function 

, „ . f d 4 k e -ik-(x-x>) 
A(x-x',M c ) = j — 



{f^C^C^x)) = {f^C^Ww) -£- 2 F, u [£( X )]F^[£(x)] = \l- 2 



5(^C„Ax)C^(x)) = (^C, 
x JJJ DxD X Dc»{d»c»{x) (fycOr) - d^(x)) + F^[£(x)](d^(x) - d v c^x))} 



expi / <fz ^c u (z)(a + Mq)c v {z) - g xm(z)j u Xm(z) c u (z) 



-9Xm{z)i„Xm(z) C v [£{z)] + XM(z)(iYdu - M) X m(z) 



(3.8) 
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For the integration over the c^-field we suggest consider the auxiliary path integral 

J(XM,XM,C V [£}) = 

Dc» {d^c u (x) (d^c u (x) - d v c„{x)) + F^[£(x)](d^(x) - cU M (x))} 



expi d z 



1 



c v {z){U + M c )c v {z)-c v {z)J v {z) 

where J v (z) is a conserving current, d v J v {z) = 0, defined by 

Jv{z) = 9 Xm{z)1uXm{z) +ju{z) 



(3.9) 



(3.10) 



and ju(z) is an external source which should be put zero finally. 

The r.h.s. of Eq.( |3.9[ ) can be represented in the form of functional derivatives with respect 
to j u (x) 



J(XM,XM,C U [£}) = 
{ d 5 ( d 







dx^ Sj u (x) \ dx^ dj u (x) dx u Sjf i (x) 



+ iF» v [£(x)] 



a 







dx^Sj u (x) dx v (x) 



rl 



x / Dd" expi / d 4 z \p;c v (z) (□ + M c )c u (z) - c u (z)J u (z) 



(3.11) 



ju=0 



At J u {z) = the path integral over is normalized to unity. The integral over c M is a Gaussian 
and the result of the integration reads 



j Dc^ expi J i 



Dc^expi / d A z ^c u (z)(a + M c )c u (z)-c u (z)J u (z) 



exp 



1 

i — 
2 



d^z dV J v (z) A(z- z', M C ) J v {z') 



(3.12) 



where the Green function A(z — z',Mc) is given by Eq. (|3.7|) . 
The functional derivatives with respect to j u (x) are equal to 







Q x fi Sj u {x) 



exp 
8 



d A z d^z' J v (z) A(z-z',M C ) J v {z') 



i I d A z- — A(x - z,M c ) J v {z) 



x exp 




dx v 5jv(x 



exp 



d 4 zd 4 z' J u {z) A{z - z', M c ) J v {z') 

' 1 f d 4 z d 4 z' J v {z) A(z - z', M c ) r(z') 



J d A z^A{x-z,M c )J^z) 



x exp 




1 

l 2 



d 4 zd 4 z' J v {z) A{z - z', M c ) J v (z') 







dx^ 8jv{x) dx^ 5j u (x) 



exp 



i- d 4 z d 4 z' J u (z) A(z - z', M c ) r(z') 



d A z D x A(x - z, Mc) Jv{z) / d A z A(x - z 1 , M c ) J u (z') 
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I d 4 z-^A(x-z,M c )Jv(z) J d 4 z'^A(x-z',M c )r(z' 



x exp 
d S i) t 



i - I I d 4 z d 4 z' J u (z) A(z - z', M c ) J u (z') 



dx^ Sj u (x) dx v Sj^(x) 



exp 



d A z d 4 z' J v {z) A(z - z', M c ) r(z') 
J d 4 z-^A(x-z,M c )J^z) J d 4 z'-^A(x-z',M c )J 
d 4 zd 4 z' J v {z) A(z - z', M c ) J u {z') 



TU (z') 



x exp 



where we have taken into account that d fM J fJi (z) = 0. 



(3.13) 



The functional J(xm,Xm,C u [£]) is then defined by 
J(XM,XM,C U [£}) = 

= \g 2 [ d A za x A{x - z,M c )[xm{z)1uXm{z)\ I d 4 z' A(x - z' , Mc)[xm(z')Yxm(z')] 



+g 2 I d A z-^-A(x-z,Mc)[xM(z) luXM (z)] J c^z' -^A(x - z' , M c ) [xm(z')Yxm(z')] 



g 2 I d A z-^A(x-z,M c )[xM(z)wi(z)] [ d 4 z' -£rA(x - z\ M c ) [xm(z')Yxm(z')] 



dx. 







+ 2gF^[£(x)} / d 4 z — A(x - z,M c )[xm{z)1uXm{z)] 



dx^ 



{-ig* l -[[ d 4 zd 4 z' [ XM {z) lvXM {z))A{z - z',M c ) [xm(z')Yxm(z')]\. 



(3.14) 



x exp 



Substituting the functional J(xm, Xm, C v [£]) in the integrand of the r.h.s. of Eq.(|3.8|) we 
express the dual-vector field condensate as a path integral over the massive magnetic monopole 
fields. 



<)(£^C^(x)C^(x) 



Z 2vr 2 




D X D X 



g 2 / d A zU x A(x - z,M c )[xm{z) 1u xm{z)\ / d 4 z' A{x - z', M c )[xm{z')Yxm(z')] 



d_ 

dx^ 
d_ 

dx h 



+</ 2 / <l ; :-^-M-r- :.3/r)[v.ul»>\.u(~)] / d 4 z' A A(x - z', M c ) [xm(z')Yxm(z')} 

_d 

dx 



1 1 d 4 z^-A(x-z,M c )[xM(z)^ X M(z)] I d 4 z' -£-A(x - z',M c )[x M (z'h u XM(z')] 



+ 2 f ,F'"-[£(.r)} I d 4 z -^-A(x - z, M C ) \xm{z)1uXm{z)] 



x expz 

-^g 2 



J d 4 z\xM{z)(i-i v d v - M) xm(z) - g Xm(z)j„Xm(z) C v [E(z)] 
\U d 4 zd 4 z< [ X M(zhuXM(z)]A(z - z',M c )[xM(z'h u XM(z')]}. 



(3.15) 



As has been shown in Refs.[3,4] due to the London limit, M 3> Mc, the non-local four-monopole 
interaction in the exponential of the r.h.s. of Eq.( 3.15| ) can be approximated by a local four- 
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monopole interactions. This reduces the r.h.s. of Eq.( |3.15| ) to the form [3,4] 

\g 2 I d A z U x A{x - z,Mc)[ X M(zhuXM(z)} [ d 4 z' A(x - z' , M c )[xm(z')Yxm(z')] 



+g 2 J d 4 z ^A(x - z, M C ) \xm(z)juXm(z)] j d 4 z> ^A(x - z', M c ) fru^'h" Xm(z')] 
-g 2 [ d 4 z^-A(x- z,M c )[ X M(zh fl XM(z)] [ d 4 z' -£-A(x - z', M c ) [xm(z')Yxm(z')] 



j dx v J dx^ 

+ 2gF^[£(x)} J d 4 z-^A(x-z,Mc)[xM(z) ll/X M(z)}\ 

x expz / d 4 z{x~ M (z)(iYdv - M)xm(z) - gx~M(zhuXM(z)C u [£(z)} 



9 2 
2M C 



[xm{z)^ u xm{z)][xm{z)Yxm{z)]Y (3.16) 



For the integration over the massive magnetic monopole fields it is convenient to decompose the 
r.h.s. of Eq.( 3.16j ) conventionally into two parts 

5^C, v (x)C^(x)) = (^C^C^ix)) 



non— string 

,2 



+ <&<W.) <?-(«)) (3.17) 



where the terms are given by 

j£-C, v (x)C^(x)) =\t~2 

. 47T Z / non-strmg Z 27T 



g 2 / d A zO x A(x - z,M c )[xM(zhuXM(z)} / d 4 z' A(x - z 1 , M c )[xM(z'h"xM(z')} 



+g 2 I d 4 z ^-A(x - z, M C ) \xm(z)7uXm(z)] j d 4 z' ^A(x - z', M c ) [xm(z')Yxm{z')} 
g 2 I d A z—A(x - z,M c )[ X m(z) Wi (z)] / efV - — A(x — z',Mc) [xm{z')^ v xm{z')) 



dx v J dx^ 



x expz j d 4 z ^XM{z){ij u d u - M)xm{z) - gXM(z)-y,sXM{z)C u [£(z)} 

9 2 - 1 

-1^2[xm{z)-1uXm{z))[xm{z)Yxm{z))^ (3.18) 



and 



/d 
d 4 z Q^ A i x ~ z, M c ) g [xm(z)juXm(z)] 



x expf j d 4 z{xM{z){iYdv ~ M)xm{z) - gXM{z)~ivXM{z)C u [£{z)] 

-[XM(z)-f u XAi(z)][XM{z)YXM(z)]y (3.19) 

13 



9 2 
2M 2 



As we would show below the non-string part of the dual-vector field condensate does not depend 
on a dual Dirac string, whereas the string part does. 

The integration over massive monopole fields is analogous to that carried out in Refs.[3— 4], 
where the r.h.s. of Eq.(3.19) has been approximated by massive monopole-loop diagrams by 
keeping only leading divergent contributions. 

In terms of the massive magnetic monopole-loop diagrams the non-string part of the dual- 
vector field condensate is determined by 



(|^C^(x)C^(x) 



d A k 

^27 tr <i ^ 



g 



non— string 
1 



16vr 2 



d 4 q 

2ttH (M 2 -q 2 -i0) 2 



1 



7rn 



TT 2 i 

d 4 k 2 



M-k + g C[S(z)] M-k + g C[£(z)] 
1 1 



n=l 



tr s 7u 7ai 

[ M M -h+gC[£(z)] M -ki+gC[S(z)] 

1 1 



jr- tr < 7 ai 



d k^i 

TT 2 i 

d 4 k nA 



:7q 



M -k 2 +gC[£(z)} M -k 2 + gC[£(z)} 
1 1 



tr < 7o 



M - k n + g 6[£{z)\ lan M-k n + g C[£(z)} 
1 1 



Ml 



M-k n+l + gC[£{z)} M -k n+l + gC[£(z)} 



16vr 2 



(3.20) 



Here we have taken into account that the contributions of the first two terms in the non-string 
part cancel each other after integration over virtual momenta. 
Integrating over k and fcj (i = 1, . . . , n + 1) we obtain [3,4] 



non— string 

oo 



16vr 2 



lf i n u 



2irH (M 2 - q 2 



2 9fMU [J 1 (M)+M 2 J 2 (M)]Y J 



n=l 



J_// 

M 2 8vr 2 



[J 1 (M) + M 2 J 2 (M)} 



32vr z 



[J^Mc) - M 2 J 2 {M C )]- 



-[J 1 (M) + M 2 J 2 (M)} 



' ^[J 1 (M)+M 2 J 2 (M)} 



M 2 8tt 2 



(3.21) 



where J\(Mc) and J 2 (Mc) are the quadratically and logarithmically divergent integrals defined 
by Eq. (|L15| ) with M replaced by M c . Since the cut-off A is much greater than Mc, we can 
replace [J\(Mc) — MlJ 2 (Mc)} by A 2 . The non-string part of the dual-vector field condensate 
is then given by 



i g 2 



|^[Ji(M)+M 2 J 2 (M)] 



non— string 



A 2 32^ 4 



1 



i g 1 



(3.22) 



M 2 8tt 2 



[J 1 (M) + M 2 J 2 (M)] 



We have represented the non-string part of the dual-vector field condensate in the form conve- 
nient for the comparison with the lattice calculations where we have set Atj = A. 
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The string part of the dual-vector field condensate is determined by the following expression 
in terms of the massive magnetic monopole-loop diagrams [3,4] 



= 9 F^[£{x)} fd 4 z/-A(x-z,M c ) 

string 7T 2 / OX^ 



9 



16vr 2 
d A k x 



d 4 k 



M — k + g C[£{z)] 



1 



+E 

n=l 

1 



a 

Ml 



1 



16vr 2 



M-h+g C[S{z)} 7m M -k+g C[£(z)} 
1 1 



M-k 2 +g C[S(z)] la2 M-k 2 + g C[£{z)} 
1 1 



r7 



M - fc n _i + <? C[£(*)] M - fc n _i + 5 C[£(z)] 



16vr 2 



TT 2 i 



M-fc n + <7C[^(z) 



Keeping only leading divergent contributions [2-4] we obtain 



/_£! 

\47T 2 



C^(x)C^(x) 



I d 4 z^-A(x-z,M c )C u [£(z)} 



string 7T Z 



Ox' 1 



8vr 2 



[Ji(M) +M 2 J 2 (M)] 



1 -4^ [Ji(M)+MV2(M)] 
Integrating by parts over z we arrive at the expression 

l-^[J 1 (M)+M 2 J 2 (M)} 

o 



d 4 z £2 (x)] A(x - z, M C ) F^[£(z)]. (3.25) 



(3.23) 



(3.24) 



Collecting all pieces we obtain the condensate of the dual-vector field 

.2 



±(£- 2 C, v {x) C^(x)) = ±£- 2 F^[£(x)] F, v [£{x)\ + 



JL 

8tt 



-,[Ji{M) + M z J 2 {M)} 



^[Ji(M) + M 2 J 2 (M)] 



M 2 8vr 2 



X I X 2 32^4 + X4 /'^S^^W]A(x-,,M c )^[^( Z )]), 



2vr 2 



(3.26) 



where the first term comes from the tree-approximation defined by the Abrikosov flux line, 
whereas the second one is fully caused by quantum field fluctuations around the Abrikosov 
flux line. Below we show that the term caused by quantum field fluctuations dominates in the 
dual-vector field condensate. 
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By using Eqs. (|1 . 19|) — ( |1 . 1 7[) and the relation G% = G/4 we can recast the r.h.s. of Eq. fl3.26Q 
into the form 

l _9^_ (xx) 

U&* c ^ x) c ^ x) ) = * & p ^ x)] p ^ x)] + "vm 

1+ M 2 M 

( M 2 3n 2 M 2 f 3« 2 \ 

X ( A*" Uto* + A* J ^^F^{8{x)]^{x-z,Mc)F^[£{z)]\. (3.27) 

The dependence of the dual-vector field condensate on the shape of a dual Dirac string enters 
through the term F^[£(x)] F^ u [£ (x)]. 

Since the confinement regime can be described by static infinitely long dual Dirac strings, we 
would make an estimate of the numerical value of the dual-vector field condensate in the approx- 
imation of the infinitely long dual Dirac string strained along the z-axis. In this approximation 
the electric field of a static dual Dirac string is given by [8] 

£(r) = e z QS^(f), (3.28) 

where r is the radius-vector in the plane perpendicular to the z-axis. The dual-vector potential 
possesses only the azimuthal component and reads [8] 

C a (r) = -^K l (M c r), (3.29) 

where K\(Mcr) is the McDonald function. The dual electric field induced by an infinitely long 
dual Dirac string amounts to 

OM 2 

E{r) = rot C(f) = - e z ^—^ K (M c r). (3.30) 

2tt 

Substituting the dual electric field Eq.( |3.30| ) in the dual-vector field condensate Eq.( |3.26| ) we 
obtain 

g 2 (xx) 

1+ M 2 M 

By using the Dirac quantization condition gQ = 2ir we bring up the r.h.s. of Eq.( |3.32| ) to the 
form 

g 2 (xx) 

1+ M 2 M 

l£ + ^^M c rK Q{ M c r)K l{ M c r)\ (3.32) 
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At r — > oo the dual-vector field condensate is defined by only the non-string part and reads 

g 2 (xx) 



U&°~u°~vh£&f±%r- (-3) 

1 + 2 q z N/ T v/ 

One can show that this value is positive. 

At r — > the dual-vector field condensate is infinite. However, as has been noted in Refs.[2— 4] 
in the MNJL model as well as in the dual Higgs model with dual Dirac strings [8-10] the minimal 
transversal distances are restricted by inequality r > 1/M a = 1/2M. Thus, at r = 1/2M we get 
the maximal value of the dual-vector field condensate, since the contributions of the dual Dirac 
string described by the terms proportional to K^{MCr) and r Kq(Mct) K\{Mcr) are positive. 

J_(ji C(< „ (l)c ,, (l) ) = M_ 

/ M 2 3g 2 M 2 Ml 3 „ , ... . . 

We should emphasize that unlike the magnetic monopole condensate the dual-vector field con- 
densate does not vanish in the close vicinity of a dual Dirac string [3,4]. 




4 Conclusion 

The evaluation of the dual-vector field condensate in the MNJL model has confirmed our state- 
ment we have pointed out in Refs. [2-4,9,10] that quantum fluctuations of the fields of monopole- 
(anti)monopole collective excitations around the Abrikosov flux lines induced by dual Dirac 
strings in the condensed phase play an important role for the undestanding of the confinement 
mechanism. We have shown that these fluctuations determine fully the dual-vector field con- 
densate independent on the shape of a dual Dirac string and non-vanishing at large distances. 

In the case of the neglect of quantum field fluctuations this part of the dual-vector field 
condensate does not appear and the condensate is defined by 



^(J^(*)C^)) = ^ (1.1) 



for the infinitely long static dual Dirac string strained along the z-axis. 

Unlike Eq.( 3.32| ) the dual-vector field condensate given by Eq. (|4.l| ) vanishes at large distances 



in the plane transversal to a dual Dirac string. This result should contradict to the gluon 
condensate having a non-vanishing part at large distances. 

When matching our expression for the dual-vector field condensate given by Eq.(|3|) with 
that calculated on lattice we would like to emphasize the absence of the perturbative term W40 
and the presence of the term proportional to M 2 /A 2 having a non-perturbative nature. The 
absence of the perturbative contribution to the dual-vector field condensate is rather clear. 
Indeed, the evaluation of the dual-vector field condensate is carried out in the non-perturbative 
(condensed) phase of the MNJL model with a non-perturbative dual-superconducting vacuum 
filled with dual Dirac strings. Thereby, the perturbative contributions cannot appear in principle. 

Then, non-perturbative contributions having the structure M^/A 4 and M^M 2 /A 4 are de- 
termined fully by the dual Dirac string. They vanish at large distances. At short distances 
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restricted from below by inequality r > 1/M a = 1/2M the contributions of these terms are 
positive. Due to this they increase the value of the non-string part of the dual-vector field con- 
densate in the vicinity of a dual Dirac string. This result is opposite to that we have obtained 
for the magnetic monopole condensate. Unlike the dual-vector field condensate the magnetic 
monopole condensate vanishes in the close vicinity of a dual Dirac string [2-4] . 

Unfortunately, we could not derive with one-to-one correspondence the contribution to the 
dual-vector field condensate analogous to that caused by a dilute-instanton vacuum in the gluon 
condensate. In the MNJL model according to the Dirac quantization this contribution should 
be of order (M^/A 4 ) 0(g & ) or (M^/A 4 ) 0(g & ). Such contributions can be found among other 
terms in the part of the dual-vector field condensate depending on the shape of a dual Dirac 
string. However, these terms vanish at large distances and do not give a dominant contribution 
to the dual-vector field condensate. Thus, we can conclude that the non-perturbative dual- 
superconducting vacuum of the MNJL model possesses partly the properties of a dilute-instanton 
gas vacuum but this part of the wave function of the non-perturbative vacuum does not play a 
dominant role. 
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